dtic  COP'- 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMBNo  0704-01-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions  searching  existing  dete  sources, 
gethering  and  maintaining  the  data  needed  and  completing  end  reviewing  the  collection  of  information  Send  comments  regerding  this  burden  estimete  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  lhe  burden  to  Department  of  Defense,  Washington  Headquarters  Services  Directorate  for  Information  Operations  end  Reports 
(0704-0188)  1215  Jefferson  Davis  Highway,  Suite  1204  Arlington  VA  22202-4302  Respondents  should  be  eware  that  notwithstanding  any  other  provision  of  lew  no  person  shell  be 
subject  to  eny  penally  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  e  currently  valid  OM8  control  number 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

25-04-2008  REPRINT 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

Low  frequency  electrostatic  waves  in  weakly  inhomogeneous 
magnetoplasma  modeled  by  Lorentzian  (kappa)  distribution 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

62 1 01 F 

6.  AUTHORS 

Bamandas  Basu 

5d  PROJECT  NUMBER 

2301 

5e.  TASK  NUMBER 

SD 

5f.  WORK  UNIT  NUMBER 

A3 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  /RVBXP 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

AFRL-RV -H  A-TR-2008- 1 132 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 
AFRL/RVBXP 

11  SPONSOR/MONITOR  S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  distribution  unlimited. 

13.  SUPPLEMENTARY  NOTES  Reprinted  from  Physics  ot  Plasmas,  Vol  15,  042108  (2008)  1)01  lit.  1063/1.2906217 

14.  ABSTRACT  l  inear  dispersion  relations  for  electrostatic  waves  in  spatially  inhomogeneous,  current-carrying  anisotropic 

plasma,  where  the  equilibrium  particle  velocity  distributions  are  modeled  by  various  Lorentzian  (kappa)  distributions  and  by  well-known  bi-Ma\we)lian  distribution, 
arc  presented.  Spatial  inhomogeneities,  assumed  to  be  weak,  include  density  gradients,  temperature  gradients,  and  gradients  (shear)  in  the  parallel  (to  the  ambient 
magnetic  field)  flow  velocities  associated  with  the  current.  In  order  to  illustrate  the  distinguishing  features  of  the  kappa  distributions,  stability  properties  of  the  low 
frequency  (lower  than  ion  cyclotron  frequency)  and  long  perpendicular  wavelength  (longer  than  ion  gyroradius)  modes  are  studied  in  detail,  and  the  results  are 
contrasted  with  those  for  the  bi-Maxwellian  distribution.  Specific  attention  is  given  to  the  drift  waves,  the  current-driven  ion-acoustic  waves  in  the  presence  of 
velocity  shear  the  velocity  shear-driven  ion-acoustic  modes,  and  the  ion  temperature-gradient  driven  modes.  Growth  rates  of  the  drift  wave  instability  and  the 
current-driven  ion-acoustic  instability  are  reduced  from  their  values  for  bi-Maxwellian  distributions  due  to  larger  ion  Landau  damping  rates  associated  with  the 
kappa  distributions.  For  the  same  reason,  excitation  conditions  for  these  two  instabilities  are  more  stringent  in  the  case  of  the  kappa  distributions  Growth  rates  of 
the  velocity  shear-driven  ion-acoustic  instability  and  the  ion  temperature-gradient  driven  instability  are  reduced  from  their  values  for  bi-Maxwellian  distribution  as  a 
consequence  of  the  reduced  adiabatic  response  of  the  electrons  to  the  electrostatic  potential  perturbation  Frequencies  of  the  drift  waves  and  the  ion-acoustic  waves 
are  also  reduced  in  kappa-distribution  plasmas  due  to  the  reduced  adiabatic  response  of  the  electrons 

15.  SUBJECT  TERMS 

Kappa  distribution 

Electrostatic  waves 

Plasma  instability 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

C.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

17 

Bamandas  Basu 

UN  CL 

UN  CL 

UNCL 

UNL 

19B.  TELEPHONE  NUMBER  (Include  area  code) 

200901 26058 


Standard  Form  298  (Rev  8/98) 

Prescnbed  by  ANSI  Sid  239  18 


DTIC  COPY 


PHYSICS  OF  PLASMAS  15.  042108  (2008)  AFRL-RV-HA-TR-2008-1 1 32 


Low  frequency  electrostatic  waves  in  weakly  inhomogeneous 
magnetoplasma  modeled  by  Lorentzian  (kappa)  distributions 

B.  Basu 

Air  Force  Research  Laboratory,  Hanscom  Air  Force  Base,  Massachusetts  01731,  USA 

(Received  24  January  2008;  accepted  17  March  2008;  published  online  25  April  2008) 

Linear  dispersion  relations  for  electrostatic  waves  in  spatially  inhomogeneous,  current-carrying 
anisotropic  plasma,  where  the  equilibrium  particle  velocity  distributions  are  modeled  by  various 
Lorentzian  (kappa)  distributions  and  by  well-known  bi-Maxwellian  distribution,  are  presented. 
Spatial  inhomogeneities,  assumed  to  be  weak,  include  density  gradients,  temperature  gradients,  and 
gradients  (shear)  in  the  parallel  (to  the  ambient  magnetic  field)  flow  velocities  associated  with  the 
current.  In  order  to  illustrate  the  distinguishing  features  of  the  kappa  distributions,  stability 
properties  of  the  low  frequency  (lower  than  ion  cyclotron  frequency)  and  long  perpendicular 
wavelength  (longer  than  ion  gyroradius)  modes  are  studied  in  detail,  and  the  results  are  contrasted 
with  those  for  the  bi-Maxwellian  distribution.  Specific  attention  is  given  to  the  drift  waves,  the 
current-driven  ion-acoustic  waves  in  the  presence  of  velocity  shear,  the  velocity  shear-driven 
ion-acoustic  modes,  and  the  ion  temperature-gradient  driven  modes.  Growth  rates  of  the  drift  wave 
instability  and  the  current-driven  ion-acoustic  instability  are  reduced  from  their  values  for 
bi-Maxwellian  distribution  due  to  larger  ion  Landau  damping  rates  associated  with  the  kappa 
distributions.  For  the  same  reason,  excitation  conditions  for  these  two  instabilities  are  more  stringent 
in  the  case  of  the  kappa  distributions.  Growth  rates  of  the  velocity  shear-dnven  ion-acoustic 
instability  and  the  ion  temperature-gradient  driven  instability  are  reduced  from  their  values  for 
bi-Maxwellian  distribution  as  a  consequence  of  the  reduced  adiabatic  response  of  the  electrons  to 
the  electrostatic  potential  perturbation.  Frequencies  of  the  drift  waves  and  the  ion-acoustic  waves  are 
also  reduced  in  kappa-distnbution  plasmas  due  to  the  reduced  adiabatic  response  of  the  electrons. 
[DOI:  10.1063/1.2906217] 


I.  INTRODUCTION 

In  collisionless  plasma,  particle  distribution  in  velocity 
space  can  depart  considerably  from  a  Maxwellian.  For  ex¬ 
ample,  in  naturally  occurring  plasmas,  such  as  plasmas  in  the 
planetary  magnetospheres  and  solar  wind  plasma,  particle 
velocity  distributions  are  observed  to  have  a  prominent  non- 
Maxwellian  (power-law)  high-energy  tail  (for  some  refer¬ 
ences  to  observations,  see  Ref.  1).  The  appropriate  distribu¬ 
tion  functions  that  can  better  model  such  particle 
distributions  are  the  generalized  Lorentzian  distributions, 
also  known  as  the  kappa  distributions.  The  kappa  distribu¬ 
tion  with  a  finite  value  of  the  spectral  index  k  has  a  power- 
law  tail  at  velocities  higher  than  the  thermal  velocity  and 
approaches  a  Maxwellian  distribution  in  the  limit  as  k— >°c. 
Typically,  space  plasmas  are  observed  to  possess  a  spectral 
index  k  in  the  range  2-6.  The  presence  of  a  substantially 
larger  number  of  suprathermal  particles,  which  distinguishes 
kappa  distribution  from  a  Maxwellian,  can  significantly 
change  the  rate  of  resonant  energy  transfer  between  particles 
and  plasma  waves.  Hence,  it  could  change  the  growth  or 
damping  rate  of  the  plasma  waves,  the  excitation  conditions 
for  instability,  as  well  as  the  rate  of  anomalous  transport 
processes  that  rely  on  resonant  wave -particle  interaction.  It 
is,  therefore,  interesting  to  study  the  stability  properties  of 
plasma  waves  when  the  equilibrium  (unperturbed)  state  of 
the  plasma  is  described  by  a  kappa,  rather  than  a  Maxwell¬ 
ian,  distribution. 

In  the  last  several  years,  plasma  waves  (electrostatic  and 


electromagnetic)  in  homogeneous,  unmagnetized  and  mag¬ 
netized  plasma  have  been  studied  by  various  authors 1,VI 4 
using  different  types  of  kappa  distributions  for  the  equilib¬ 
rium  state.  In  this  paper,  we  concern  ourselves  with  the  elec¬ 
trostatic  waves  in  spatially  inhomogeneous,  current-carrying 
anisotropic  plasma,  where  the  equilibrium  particle  velocity 
distributions  are  modeled  by  various  kappa  distributions. 
Spatial  inhomogeneities,  assumed  to  be  weak,  include  den¬ 
sity  gradients,  temperature  gradients,  and  gradients  (shear) 
in  the  parallel  (to  the  ambient  magnetic  field)  flow  velocities 
associated  with  the  current.  Such  equilibrium  plasma 
configuration  is  representative  of  many  space  and  laboratory 
plasmas.  We  consider  three  specific  forms  of  the  kappa 
distributions,  namely,  kappa-Maxwellian,  product 
bi-Lorentzian,1  and  bi- Lorentzian.1  We  first  present  the  full 
dispersion  relations  for  electrostatic  waves,  then  present  the 
versions  of  the  dispersion  relations  that  are  more  suitable  for 
the  study  of  low  frequency  (lower  than  ion  cyclotron  fre¬ 
quency)  waves,  and  finally  concentrate  on  the  analysis  of  the 
low  frequency  and  long  perpendicular  wavelength  (longer 
than  ion  gyroradius)  modes.  In  particular,  the  stability  prop¬ 
erties  of  drift  waves,  current-driven  ion-acoustic  waves  in  the 
presence  of  velocity  shear,  velocity  shear-driven  ion-acoustic 
modes,  and  ion  temperature  gradient  driven  modes  are  ana¬ 
lyzed  in  detail.  Such  plasma  modes  are  commonly  observed 
in  inhomogeneous  plasma  and  are  responsible  for  anomalous 
effects,  such  as  diffusion,  thermal  conduction  and  resistivity. 
We  also  include  in  our  presentation  the  corresponding  dis- 
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persion  relations  for  the  bi-Maxwellian  distribution  and  their 
analysis.  The  reasons  for  including  the  previously  known 
results  for  the  bi-Maxwellian  distribution  are:  First,  to  check 
the  correctness  of  the  results  for  the  kappa  distributions  by 
using  the  fact  that  the  kappa  distribution  goes  over  to  the 
Maxwellian  distribution  in  the  limit  as  second,  to 

compare  and  to  illustrate  the  distinguishing  features  of  the 
kappa  distributions. 

The  paper  is  organized  in  the  following  way:  In  Sec.  II, 
we  describe  the  general  mathematical  formalism  for  the  deri¬ 
vation  of  dispersion  relation  in  weakly  inhomogeneous  mag¬ 
netoplasma.  In  Sec.  Ill,  we  present  the  dispersion  relations 
for  the  different  equilibrium  distribution  functions  described 
in  Sec.  II.  Section  IV  is  devoted  to  the  analysis  of  the  low 
frequency,  long  perpendicular  wavelength  modes  for  the  dif¬ 
ferent  equilibrium  distribution  functions.  The  results  are 
summarized  in  Sec.  V. 


II.  GENERAL  MATHEMATICAL  FORMALISM 


A.  Equilibrium  state  (spatial  inhomogeneity 
along  the  x-direction) 


The  dynamics  of  nonrelativistic,  collisionless  plasma  is 
determined  by  the  Vlasov  equation 


S  _  qa  I  _  v  X  B 

— +  v  V+  —  E  + - , 

St  mn\  c  J  S\ 


/„(r,v,r)  =  0,  (1) 


where  fa  is  the  single-particle  distribution  function,  qa  is  the 
charge,  and  ma  is  the  mass  of  the  plasma  constituent  a  (a 
-e  for  the  electrons  and  a-i  for  the  ions),  while  E  and  B  are 
the  electric  and  the  magnetic  field,  respectively.  The  distri¬ 
bution  function  for  the  steady  equilibrium  state  of  spa¬ 
tially  inhomogeneous  plasma  immersed  in  a  uniform  mag¬ 
netic  field  B0  obeys  the  time- independent  Vlasov  equation 

|v.V  +  —  vXB0-  7- j/aoO-.v)  =  0  (2) 

\  rrio/c  d\) 


and  it  can  be  constructed  from  the  constants  of  motion  of  the 
charged  particles.  If  we  adopt  a  cylindrical  coordinate  system 
in  velocity  space  with  its  z-axis  parallel  to  B0,  so  that 
vx=v1cos  (p,  vy=v1  sin  (p  and  vz  =  v-B0/Z?0=i; ,  then  the 
constants  of  motion  are  jc +vy/Qcr  y-vx/fla,  + 

and  v  .  Here,  is  the  cyclotron  frequency  of 

the  charged  particle  species  a.  We  assume  that  the  plasma  is 
inhomogeneous  only  along  the  jc  direction.  Then,  the  most 
general  is  a  function  of  fsx+Uy/fl*,  l>^,  and  That 
/aO=/ao(f»ui»i;ii)  is  a  solution  of  Eq.  (2)  can  be  verified  by 
direct  substitution.  If  the  spatial  gradients  are  weak, 
fao(€>v\  ,U||)  can  be  expanded  in  a  Taylor  series  about  £=0. 
Thus,  retaining  only  the  terms  that  are  linear  in  the  gradients, 
we  have 

/ao(£wi.t'n)  s  [1  +  &L(v\,vy\FM(,v\,v  ),  (3) 

where  /rao(«i.«ii)*/ao(^=°»i'i»i'ii)  and  1  IL{v\,v^ 
which  contains  the  essential  features  of 
the  inhomogeneous  plasma  including  density  gradient,  tem¬ 
perature  gradient,  and  shear  in  the  flow  velocity.  In  the  fol¬ 
lowing,  we  consider  the  different  specific  forms  of 


mentioned  in  the  Introduction,  and  derive  the 
respective  expressions  for  1/L. 

(a)  Bi-Maxwellian  (BM) 


”a0(j) 


Xexp 


pj  4(g) 


(4) 


where  ««>(£)  =  v  and  Va0(^)  is  the  inhomogeneous 

parallel  flow  velocity.  Then,  according  to  Eq.  (3),  we  find 


L~  Ln 


1 

LoaW 


-0a 


1 


ei 


-I  1  Ul  ^ 


-Hn  Oa 


(5) 


and 


”«() 


(6) 


Here  \/Lna=(d\nnM/d0^o,  MLm=(d\n  6a /^)^0, 
\IL6a  ={d  ln  Oal/d(;)f=0,  and  Va0=(dVM/dQ(m0-  It  is  to  be 
understood  that  /i„o,  VM,  0 and  0al  appearing  in  Eqs.  (5) 
and  (6)  denote  their  values  at  £=0.  We  have  suppressed  the 
arguments  for  simplicity  in  notations.  The  function 
>P||)  >s  normalized  to  the  density  nM,  and  the  thermal 
speeds  0&  and  0a^  are  related  to  the  particle  temperatures, 
Ta  and  TaX,  by  &M=2TMlma,  &a  =2TaJm„,  where  the 
definitions  of  T-oii  and  Tnl  are 

>'aOTai  =  2mna  J  dv1dv,,v1(v  -  V^F^v^v  ), 


nrtTa±  =  ™na  J  dv xdvv[)\F . 

(b)  Kappa-Maxwellian  (KM) 

n«o(£)  r(#f) 


(7) 


^vU)  = 


WaAMJd  Kmr(K-  1/2) 


X 


- IP 

R 

^  r-j 

_ l 

-K 

1 

v\ 

*<£<(£)  J 

exp 

L"«L<eJ 

(8) 


where  T  is  the  gamma  function  and  =u  -  Va0(^). 

Then,  according  to  Eq.  (3),  we  find 


1  1  1 

L  Lna  L#a\ 


+  2 


1  uli\  u«V'ao\l  ,  “1 

t  et  k*. 


-1 


(9) 


and 
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pKM(  2  s_  "oO _ _ 

°°(  1  v312!?  6MKmY{K- M2) 


Ol  UCt\\  ' 


X  1  + 


4a\\ 


exp  - 


(d)  Bi-Lorentzian  (BL) 

HgoU) 


(10) 


/Sotfcwl.Wy)  - 


ru) 


As  before,  ««<),  V^,  0^,  and  0al  in  Eqs.  (9)  and  (10)  refer 
to  their  values  at  £=0.  The  function  F^(i;x  ,tf|j)  is  normal¬ 
ized  to  the  particle  density  n^,  while  0 ^  and  da x  are  related 
to  the  particle  parallel  and  perpendicular  temperatures, 


^/20L(a<ua*1,2r(*-i/2) 

u2  ujiQ 

k(L(&  k&M). 


(14) 


where  w^f)  =  v\  ~  F<*o(f)-  Then,  according  to  Eq.  (3),  we 
have 


01a^  =  2TaXl ma  for  k>3/2 

(c)  Product  Bi-Lorentzian  (PBL) 


'  o0 


l  + 


l  + 


"2J& 


«Ax(Q 


\k-3)I  K](TMlma)  and 

1  1 

1 

2  2(i 

1 

L  Lnot  ^Sor 

+ 

x( 

1  u\ 

1 

r(*„ 

+  D 

l 

L#ax  0^x 

T  ^ 

Lm  0^ H 

'ii+  1/2) 

x( 

v2 

1  +  1 

2  \  -1 

+^f) 

-<*!+!) 

4  a  vai) 


(15) 


-(*i+l) 


(11) 


where  =  vl{- Vao(£).  Then,  according  to  Eq.  (3),  we 

find 


and 


«o0 


rbL/  l  \  _ 

n  .A/2 


r  M 


*1/2ru- 1/2) 


,  r2 

X  1  + 


V 


><#■„ 


-<*+l) 


(16) 


I__L 

L  ~  U, 


1 


2  2(kl  +  1)  1  u2 

+ 

-'Sal.  K 


X  1  + 


*  Or 


x  l  +  ■ 


-l 


Lfial  O^a  1 

2(*||  +  1 )  /  1  “all  +  ua  YaO 

Lsa\\  0^ ||  0^, 


02) 


and 


r(«  *i) 


+ 1/2) 


X  1  + 


“(*H+1) 


1  +• 


*  01 


-<*.+!) 


03) 


As  before,  n a0,  V^,  0a,„  and  6al  in  Eqs.  (12)  and  (13)  refer 
to  their  values  at  £=0.  The  function  F^L(y^  ,v  )  is  normal¬ 
ized  to  the  particle  density  n a0,  while  ^arll  and  0a^  are  related 
to  the  particle  parallel  and  perpendicular  temperatures, 
defined  by  Eq.  (7),  as  0^,  =[(2k  - 1  )f  k^T^I ma)  and 
0^^  =  [(/c1-l)//c^](27’crX/ma),  for  k>  1/2  and  #Cj_ >  1 . 
Here  we  have  allowed  the  possibility  of  different  values  of 
the  spectral  index  in  the  parallel  and  the  perpendicular  direc¬ 
tions.  If  k  =k±  =  k,  (13)  corresponds  to  the  product  bi- 
Lorentzian  listed  by  Summers  and  Thome. 


where  n^,  V ^  8^,  and  0aX  in  Eqs.  (15)  and  (16)  denote  the 
values  of  these  quantities  at  £=0.  The  function  ,t>n)  is 

normalized  to  particle  density  n^,  while  0a  and  8al  are 
related  to  the  particle  parallel  and  perpendicular  tempera¬ 
tures,  defined  by  Eq.  (7),  as  0* ,  =[(2/c-3)//c](7^  Im a)  and 
0i±=[(2K-3)lK](Tal/ma),  for  k>3/2. 

It  may  be  venfied  that  for  all  the  choices  of  ,u(!) 

the  equilibrium  values  of  the  density  and  the  temperatures 
are:  na(x)=na0(\ +xl  Lna),  Tall(x)  =  7^,(1  +*/L7a!),  and 
7a^(jf)  =  TaX(  1  +.r/ L7aX)  to  the  lowest  order  in  ( x/L ),  where 
LTal  =LSa\\/2  and  LTa±-L0aLl2  are  the  scale  lengths  of  the 
temperature  gradients.  Charge  neutrality  of  the  equilibrium 
state  requires  Ln€—Lnv  There  is  an  equilibrium  current  den¬ 
sity,  which  is  a  combination  of  the  diamagnetic  current  den¬ 
sity  and  the  current  density  due  to  flow  velocity  along  the 
ambient  magnetic  field  B0.  The  self-consistent  magnetic  field 
due  to  the  equilibrium  current  density  is  assumed  to  be  neg¬ 
ligibly  small  compared  to  the  main  ambient  magnetic  field 
B0. 

For  simplicity  of  notations,  we  assumed  above  that  the 
spectral  index  k  has  the  same  value  for  both  electron  and  ion 
distributions.  However,  the  analysis  in  the  following  sections 
can  be  easily  generalized  to  allow  different  values  of  k 
for  the  two  charged  populations.  In  Fig.  1  we  have  presented 
G(u)  =  (2'rr)3/2(VTail/naQ)jdv1v1Fa()  as  a  function  of 
u  =  (v  -VM)/VTai],  where  VTal=(TMlma)vz,  for  the  four 
model  velocity  distributions  described  above.  It  may  be  veri¬ 
fied  that  G(u)  for  bi-Lorentzian  is  same  as  that  for  kappa- 
Maxwellian.  As  is  well  known  in  plasma  physics,  the  func¬ 
tion  G(u)  plays  important  role  in  determining  the  stability 
properties  of  plasma  waves. 
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G(u) 


FIG.  1.  Comparison  of  G{u)  =  (2'rr)i/2(VTa{[/na0)  idvLvLFM  vs 
u  m(vt-Vc^)/VTaa.  Solid  curve  is  for  bi-Maxwellian  distribution,  dashed 
curve  is  for  kappa-Maxwellian  (k=3)  distnbution,  and  dashed-dotted  curve 
is  for  product  bi-Lorentzian  (^=3)  distribution.  G{u)  for  bi-Lorentzian  is 
the  same  as  that  for  kappa-Maxwellian. 


B.  Perturbed  state  (electrostatic  perturbation) 
and  dispersion  relation 

We  now  consider  electrostatic  perturbation  so  that 
the  perturbed  electric  field  £i(r,/)  is  given  by  £i(r,r) 
=  -V</>1(r,r),  where  <£i(r,/)  is  the  perturbed  potential.  The 
equations  to  be  solved  are  the  linearized  Vlasov  equation 
for  the  perturbed  distribution  function  /al(r,v,/)  and  the 
Poisson  equation  for  ^(r.f).  They  are 


1  -X S  f d-  ^ - 

a  *=-*  J  OJ-  k  U|j  +  /2ila 

( k,  1  d  d  \ 

xljtl+*^“2"nwr"0<1’1-t')=0- <l9) 

where  7„(/0  is  the  Bessel  function  of  the  first  kind  and 

In  this  paper  we  wish  to  study  in  detail  the  low  fre¬ 
quency  waves,  (coa7k  VY^)  1  in  inhomogeneous 
plasma,  where  coa=  co-k^V^.  The  reduced  dispersion  rela¬ 
tion,  which  is  adequate  for  the  study  of  such  low  frequency 
waves,  is  obtained  by  first  rewriting  Eq.  (19)  as 


=  0  (20) 


after  using  nQa=co-k  v  +nfla-((o-k  u  )  and  the  Bessel 
function  identity  £„./“=  1,  and  then  keeping  only  the  n- 0 
term  in  the  summation.  The  result  is 


1  - 


2[1  -4(fia)]jT- 

ovx 


a)  —  A'||U|| 


\naL  dvj 


Fm{v\,v,)  =  0. 


(21) 


—  +  v  •  V  +  (v  X  B0)  • 
dt  rric/c  d\ 


fa\(r,\,t) 


-  —  V  4>\  •  ~fa o(r,v)  =  0 


and 


It  can  be  venfied  from  the  results  presented  in  Secs.  Ill  and 
IV  that,  under  the  conditions  ( coa ,  *  VrJ<n„,  the  n¥=0 
terms  in  the  summation  in  Eq.  (20)  are  negligibly  small  com¬ 
pared  to  the  /t  =  0  term,  for  all  the  equilibrium  velocity  dis¬ 
tributions. 

In  the  next  section,  we  present  the  dispersion  relations 
that  are  obtained  for  the  different  equilibrium  distribution 
functions  after  substituting  the  specific  expressions  for  l/£ 
and  F ^  into  Eqs.  (19)  and  (21)  and  evaluating  the  velocity 
integrals. 


V20,(r,/)  =  - 


(18)  III.  DISPERSION  RELATIONS  FOR  THE  MODEL 
DISTRIBUTION  FUNCTIONS 


In  solving  Eqs.  (17)  and  (18)  we  adopt  the  ‘local 
approximation”  **  and  assume  that  the  perturbed  quantities 
have  space-time  dependence  of  the  form,  Aj(r,/) 
=  Aj(k, (o)exp(ik-r- icot),  where  k  =  (0 ,k±yk).  The  “local 
approximation”  requires  k±  ^>dldx^>  1  !L  and  pJL<^  1, 
where  pa  is  the  gyroradius  of  the  particle  a  and  L  represents 
the  typical  scale  length  of  spatial  inhomogeneity.  As  is  well 
known,  the  “local  approximation”  retains  the  leading-order 
effects  of  the  spatial  gradients.  Using  Eq.  (3)  for  /^(^v) 
and  solving  Eqs.  (17)  and  (18)  by  the  standard  procedure,1 
the  linear  dispersion  relation  for  electrostatic  waves  in  inho¬ 
mogeneous  plasma,  under  the  “local  approximation,”  is  ob¬ 
tained  as 


(a)  Bi-Maxwellian  (BM):  When  F^v2  %v)  is  a  bi- 
Maxwellian,  given  by  Eq.  (6),  and  l/£  is  given  by  Eq.  (5), 
Eq.  (19)  yields  the  dispersion  relation 


r„<« 


k2(f  “ 

*  ua\\n=-*> 


1 

20, 


a  N  ^na 


■'Oct 


+ 


r. 


zBM(^„) 


^  +  ”0 


^f’of'Ootw 


dfn 


k  ea 


zBM(0)  i 


i 

2  \  *,  O 

=  0,  (22) 


where  jt)2a=47r<7^ia0/m„,  fia=k\ ^i/(2fi2),  co^co-kV^, 
£na=(toa+ntoa)f(k\\&di\}’  and  ZBM  the  well-known  plasma 
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dispersion  function  associated  with  the  bi-Maxwellian  paral¬ 
lel  velocity  distribution,  defined  by 


ZBM(s) 


-f  r 

yj7TJ  -  oc 


ds 


exp(-  s2) 


Im  s  >  0 


s-s 


(23) 


and  by  the  analytic  continuation  of  Eq.  (23)  for  Ims^O. 
Using  &a\r2TJm<r  &L±=2T<*±/ma'  Lm=2Lt^  and 
L8al  =  2LTal ,  where  \/LTa{\=(d  In  TJdg)^  and  1  ILTal 
=  (d  In  Ta  Jdg)^  we  rewrite  Eq.  (22)  in  the  more  familiar 
form, 


i-2 


k2\2 

K  Ar>, 


2  \\{ba) 


1  r; 

1  ~~Vcr  +baVc,±-fT 


+  —nCta 


z™(L 


v2  kV, 


II  v  Tail 


+ 

2 


0 

k  n„ 


-  Vdi' 


w*„(w0  +  ;  iO„) 
2*^, 


— ZBM(£no)  f  =0.  (24) 


Here,  ^=o)a/(^&„Vral). 

(b)  Kappa-Max wellian  (KM):  When  /7a0(r‘  , ^u)  is  a 
kappa-Maxwellian,  given  by  Eq.  (10),  and  1/L  is  given  by 
Eq.  (9),  Eq.  (19)  yields  the  dispersion  relation 


i-2^2  r„03J 

a  *  ^oil«=-0C 


2H„ 


*  ^  O 


1 


VJOT 

.KM/ 


_L  + 

k-eoii 


Lea.  F 


If, 

•«  2\  k  a. 


it  0. 


aja  +  /i(le 


1=0. 


(27) 


where  Pa=k\<?aJ(2tfa),  ^a=Uoa+n(lJ/(k  BJ.  and  Z™ 
is  the  plasma  dispersion  function  for  the  kappa-Maxwellian 
distribution,  defined  by 


Here  7(477(7^/1^)),  7j0^'—Lna/Lj0f.  Vai.  k.^1 

u)t,a-k1T^I{maCtaLna),  („a=  (<5a+n(la)l  ( \2k^Vr<^), 

^niba)=^niba)exPi~ba)’  where  /„  is  the  modified  Bessel 
function,  and  T'  is  the  derivative  of  Tn  with  respect  to  its 
argument  ba=kr1Tal/ (mj.\2n). 

The  reduced  dispersion  relation  for  the  low  frequency 
waves  in  bi-Maxwellian  plasma,  derived  from  Eq.  (21),  is 


2(o‘ 


2co : 


i  +  2  f^[i  -  r  „(&,)]  -  2  t rfWa) 


k  dt 


k  Or 


all 


kJLi  1  1  2/3aF'\zBM(U 


20. 


LdaL^Q'  h&a 


+  i 

2  \  ku(la 


(25) 


or,  in  terms  of  ba,  ,  r/a  ,  and  co*a,  is 


i-i’o(U  0(O 


,  v  1  *  o yua>  v1 

1  +  Zj  ~ 7771 ~  2j 


kz\ 


Da  1 


2x2 


kz\ 


■Doll 


1  r0 

1  -~Vai  +ba7]a  — 

2  ro  /  \  2k  VTai\ 


Z™(U 


^2^V^J^ZBM(U 


=  0. 

(26) 


Z™(s)  = 


ru) 


7rl/2A-,/2r(/f-  1/2) 


X 


(5  -<?)(!  +  s2/*)*’ 


Im  s  >  0 


(28) 


and  by  the  analytic  continuation  of  Eq.  (28)  for  lm  It 
was  first  introduced  and  studied  by  Hellberg  and  Mace."  Us- 
mg  &tt=[(2K- ?>)! K\{TMlma),  &La^2TaJma,  LBa=2LTa{{, 
and  Ltial=2LTaX,  we  rewrite  Eq.  (27)  in  terms  of  ba , 

?7ai,  and  w*a  as 


i-2 


i 


k2kn„  2k-  3 „=_, 


2k  ^  ,  J  /  2/r-  3 ' 1,2 

2  Wl 


2/c 


.  1  xK  F«i  ^  T*  n 

1  -“7af  +  OaJ7a±-prjw*a+— «Oa 


x  ^(Q  f  1 

ak  VT->,  2 


k  V' 


1  -  — - 


rTal‘ 

<*>*a(<*>a  +  nila) 


-  Vail' 


2  k2V2 


k  nc 

d 

od<a 


zr«jf=0  (29) 


where  now  s*a=[2*/(2*-3)]1/2(wa+/7n a)f(\  2k  VT&). 

The  reduced  dispersion  relation  for  the  low  frequency 
waves  in  kappa-Maxwellian  plasma,  derived  from  Eq.  (21), 
is 
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i  +  X  ^[i  -  r0(/y]  -  X  |%r0(/g 

n  *  n  K  tflj  | 
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2(1  „ 
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1  +  n\z™K) 


“ 6a\\ 


L#al  r0 


hd„ 


k  v' 


+  T  — 


oO 


*1^1  ^  \  d  ^KM/  k\ 

£>2  /  J  *  * 


k  fla  (\J^Qa\\k  &“a\  f  d^a 


or,  in  terms  of  ba, ,  77^,  ?7al,  and  a)*0,  is 


=  0 


(30) 


J  +  ^  1  ~~  r 0 (ba)  r o(^a)  2  AC 


Z-2\  2 

*  ADal 


A'2\L  2  AC -3 


2ac  -  3 
2k 


1/2 


/  1  r'\ 

X  1  +  IT  2^(0 

\  2  l0/ 


+V^. 

2 \  *„  n* 


*7«! 


2k2  V2 

LK  VTa\\ 


— ZKM(c*) 


=0. 


(31) 


Here,  s^=[2/c/(2/c— 3)]1  2o>al  (y2kltVTo, ).  In  the  limit  as 
*-».  <*-£„«.  sS-*fcr  Z-M-ZBM,  and  thus  Eqs.  (29) 
and  (31)  become  identical  to  their  bi-Maxwellian  counter¬ 
parts  Eqs.  (24)  and  (26),  respectively. 

(c)  Product  Bi-Lorentzian  (PBL):  When  F^i v\,v  )  is  a 
product  bi-Lorentzian,  given  by  Eq.  (13),  and  ML  is  given 
by  Eq.  (12),  we  use  the  series  expansion  of  J2  given  by18 


.2/  .  Z  (-ir|n  r(2p  + 1)  /  pX* 

*  ^|(p-W)!(p!)Ir(p  +  |»|  +  l)\2) 

in  Eq.  (19)  and,  after  carrying  out  the  integrations  in 
-v±  space,  obtain  the  dispersion  relation 


2  +x  +» 

'-2jj£2  2 

a  K  ua\  n=-^  p=}n\ 


«J,j(KL-p)(k21e1nj 


k±&^i  w0,+  nll0\  d  nBL,  K 
(q" 


jpfp 

K  rt& 


=  0 


(33) 


for  K±>p.  Here  sj0=(5a+nfla)/(*||dot|). 

(-i)HM  rC2p+o 

SH*  p!(p-|i.D!r(p+H  +  i)’ 

and  ZPBl  is  the  plasma  dispersion  function  for  the  product 
bi-Lorentzian  distribution,  defined  by 

r(Ac„+i)  r  ds 


Z^BL(q)  = 


\77/f1/2r(/f1+i/2) 
Im  q  >  0 


f 


(j-^d+f2/^1, 

(35) 


and  by  the  analytic  continuation  of  Eq.  (35)  for  Ims^O. 
Apart  from  the  multiplication  factor,  this  function  is  the  same 
as  the  dispersion  function  introduced  and  discussed  by  Sum¬ 
mers  and  Thome.1  Using  0^=[(2aC||- O/ACjj^r^/mJ,  0lal 
=[(*±-1)/ K±](2Tal/mJ,  L8al=2LTM,  and  LAnL  =  2LTaL 
we  rewrite  Eq.  (33)  in  terms  of  ba,  77^,,  rjal,  and  a)*a  as 


Y  1  2 AC;|  y  Y  J  (*1  P)  l  K 

^  .?a2  91^  "  8\n\'P 

2*  1  n=-oc  ^=^1 


(36) 


where  now  S^=[2/c 7(2*,,-  \)]v2(a>a+nSlJ/(\!2kaVTal). 

The  reduced  dispersion  relation  for  the  low  frequency  waves  in  product  bi-Lorentzian  plasma,  derived  from  Eq.  (21),  is 


1  _  V  I^V  //  1 1  (  ac  -  p  +  1 )  /  k\  6ra^  V  y  -wna y*  ^  ^(ac  -  p)  I  k16ral  V 
n*  )  Ufi2/  „k2<?M£>8o'p  r (kx)  Uo2/ 

Ak_JL(  1  1  ,  2p  \z^K)  1/  kL  vM  kX  <  \  d  ZPRI)  ] 

1  2Cla\Lna  Lenll  Lgal )  kn0M  2  \  k  no  d*Ka  *  al 


(37) 


or,  in  terms  of  ba ,  77^,  7jal,  and  o)*a,  is 
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1  +0C 

“  2  7 2\  2  ^  8 o,/ 
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"a  )  2j  ,2.  2 


(ki-DHkx)  V  2k 
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a  1  (  L  V*  ^*a^a  \  /V 


=  o 


(38) 


for  k±>/>.  Here  s*= [2 *||/ ( 2*ii - 1  )]v2ioa/( \^.ktVTM).  The 
condition  K^>p,  which  arises  from  the  integral 

does  not  present  any 
difficulty  in  using  Eqs.  (36)  and  (38)  for  the  study  of  long 
perpendicular  wavelength  (b<*<^Cl)  modes,  including 
finite-ba  corrections,  if  k±^2. 

At  this  point,  let  us  verify  that  Eq.  (36)  — >  Eq.  (24)  and 
Eq.  (38)  — *  Eq.  (26),  in  the  limit  as  /cN,  #cx  — »►  o°.  To  do  this, 
we  note  that  ZBBL-ZBM  in  the  limit  as  k j|— >°°, 

and  that  Kpir(K±-p)/V(K1)  — ►  1  in  the  limit  as  >oc.  As  a 
result,  Eq.  (36)  becomes 
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(39) 


we  once  again  use  the  series  expansion  for  j\,  given  by  Eq. 
(32),  in  Eq.  (19)  and  after  carrying  out  the  integrations  in 
v  -v±  space  obtain  the  dispersion  relation 


0  +<»  +oc 


i  _  y 2oipa  y  yg  _  Y 

VM 


4ft: 


2  n 


i 


a  ' 


1  2  p 

; —  +  - - 

'6a  1 


1 


$a  1  1  K  uct 


z  l(qk  ) 

h  a 

K  C7- 


<^i  ( K~P\ «ft« 


17  ail 


^BL 

l  0  K'p 
K  OcA\ 


k  v: 


C.(C)  +  :'-7 


rrO 


a-  n„ 
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(42) 


for  /c >/?.  Here  s^a=(iu«+/ina)/(/:|  0tfi)  and  we  have  intro¬ 
duced  the  generalized  dispersion  function  Z^,  with 
(7=0, 1 ,2,3,.. .,  for  the  bi-Lorentzian  distribution,  defined  by 


7blm  _  »  ru)  r  ds 

*■’  N77/c1/2r(/c-l/2)J.x  (*-?)(!  H-*2/*)*-*’ 


Next,  multiplying  Eq.  (32)  by  y^exp[-mai/^/(27al)]  and 
integrating  from  i/1=0  to  i/1  =  00,  we  find 


X  «|«|,p(  ~f  )  =  In(ba)exp(-  ba)  =  r n(ba)  (40) 

H«l  V  2  7 


Im  <;  >  0 


(43) 


and  by  its  analytic  continuation  for  Ims^O.  We  note  that 
o  is  identical  to  the  dispersion  function  7™  defined  by 
Eq.  (28).  This  has  to  do  with  the  fact  that  fdv±v^F™ 
=fdv±v1F^).  Differentiating  Eq.  (43)  with  respect  to  s  and 
then  performing  integration  by  parts  we  find  the  relation 


and,  hence, 
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(44) 


It  is  then  evident  that  Eq.  (39)  is  identical  to  Eq.  (24)  by 
virtue  of  Eqs.  (40)  and  (41).  In  a  similar  manner,  it  can  be 
verified  that  Eq.  (38)  — +  Eq.  (26),  in  the  limit  as  K\yK±-+™. 

(d)  Bi-Lorentzian  (BL):  When  Fa0(i/^,U||)  is  a  bi- 
Lorentzian,  given  by  Eq.  (16),  and  1/L  is  given  by  Eq.  (15), 


which  can  be  used  to  generate  Z^  for  <7=  1,2, 3,...  from 
Z^=0.  As  in  the  previous  cases,  we  use  0^,  =  [(2k-3)/k] 
X(rdl/ma),  &al=[(2K-3)/ K](Tallma)y  Lm=2LTw,  and 
L6a±  =  2LTa±  to  express  Eq.  (42)  in  terms  of  7}al ,,  rjnl , 

and  We  obtain 
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The  reduced  dispersion  relation  for  the  low  frequency  waves  in  bi-Lorentzian  plasma,  derived  from  Eq.  (21),  is 
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(45) 
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Here  <;*=[2k/(2k-3)]1/2w0/(\2A:  Vroil).  It  can  be  verified  by 
following  the  analysis  given  for  the  product  hi-Lorent/ian 
that  Eq.  (45)  — >  Eq.  (24)  and  Eq.  (47)  — >  Eq.  (26)  in  the  limit 
as  k— *oo. 


Zbm(£,)  =  Nirexpf-  £2)  -  ~(  1  +  ^2  +  >  1 

(49) 


in  Eq.  (26)  and  then  keeping  the  leading  terms,  we  find 


IV.  LOW  FREQUENCY  AND  LONG  PERPENDICULAR 
WAVELENGTH  1)  MODES 

In  this  section,  we  present  analysis  of  the  dispersion  re¬ 
lations  for  the  low  frequency  waves,  in  order  to  illustrate  the 
differences  between  the  equilibrium  velocity  distributions. 
For  analytical  tractability  we  restrict  ourselves  to  the  situa¬ 
tion  where  |<o,/(\2 k  Vt>  )|  <C  1  and  \a)j/(\2k  Vr(11)|  1.  Fur¬ 

thermore,  we  consider  long  perpendicular  wavelength  modes 
such  that  be** 0  and  fc;<l.  These  are  the  limiting  conditions 
under  which  drift  waves,  ion-acoustic  modes,  and  ion  tem¬ 
perature  gradient  driven  modes  are  excited.  We  start  this  sec¬ 
tion  with  the  review  of  the  previously  known  results  for  the 
bi-Maxwellian  case  and  then  discuss  the  various  kappa- 
distribution  cases. 

(a)  Bi-Maxwellian  (BM):  We  refer  to  Eq.  (26).  Using 
r0(*,)~l,  r0(bt)st\-bh  k\Df u<C  1  (quasineutrality  condi¬ 
tion),  and  the  expansions  8 


ZBM(£)  =  iViKl-e+  --)-2fc  1--C  + 


(48) 


k  v'n  „  .,  “*<• 

1  +—  TT+klp;-—  [1  -bj(  1  +  T)iL)] 

k,\  li,  w,- 


k  V'n  <X>*' 

1-“-—  U  +  tt  ~b,Vu) 

k  llj  (X): 


k;c; 


+  i\  - 


2  knV 


Te  L 


A.'n  ilf 


1  -  -T) 
2 


+  i\  — 


2\kV- 


b,r), 


Ti 


k2c ? 


k  V’ 


1-— -f -  — 
k  a  at, 


i  -- 


V, 


—  +  ~Vi—  exp(-— J-T- 

_j  2  o)j  J  \  2k  Vm 


_*> 


=  o. 


(50) 


where  p^=c^/Ylj,  c]=Tfll/mh  and  we  have  used 

We  first  consider  the  drift  waves  (w,  —  ^,).  which  are 
realized  when  is  so  small  that  k^c^o)^.  Drift  waves  have 
been  discussed  extensively  in  the  plasma-physics  literature. 
Here  we  discuss  some  of  the  salient  properties.  The  relevant 
dispersion  relation  that  follows  from  Eq.  (50)  is 
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k  V' 

l+“f 

k  iL  a ), 


2  kV' 


k  VL 


eO 


Tel  L 
-2 

OK  60 
77  i 


K  iL  Q),  \  l 


+ 1  \  o  w 


8'"W; 


s? 


“p(-i^vr)=° 


riii 


fV  vm 

and  its  approximate  solution  is 

Re  l  -6,(1  +  Jfcjjw*,, 


(51) 


(52) 


/  77  .  Re  60.  I 
Imw.s  Xj-A(ki>kl)— U 

Ki\vTet\  { 


'-~Ve  H* 


-  W*, 


^Te  J  ^ / 


Xexp 


^r/ii  V  \  IfcuV'T'/ii 
Re  60; 


V2 k  V- 


7/11 


where 

*(*±.*n)  = 


2  J1  * 


i  +(My(v^n,)+*ip; 


(53) 


(54) 


We  recall  here  that  60 i=oj-k  V®  and  aie  =  a)-  k^V^  aj{-  kl{Vd, 
where  Vd =  Vtf0-  V,-0  is  the  drift  speed  of  the  electrons  relative 
to  the  ions.  The  solution  describes  drift  waves  (in  a  reference 
frame  moving  with  the  ion  flow  velocity)  modified  by  the 
presence  of  temperature  gradients  (17*11,17/1)  and  sheared 
electron  flow  velocity  (V'0).  Typically,  \(k ±/ k)(V’eQ/ fie)\  is 
very  small  compared  to  unity  and  may  be  neglected.  The  first 
term  in  Eq.  (53)  represents  resonant  interaction  of  electrons 
with  the  waves  (Landau  resonance)  and  it  gives  rise  to  insta¬ 
bility  when  the  proper  condition  is  met.  The  last  term  in  Eq. 
(53)  represents  resonant  interaction  of  ions  with  the  waves 
and  leads  to  Landau  damping  or  growth  depending  on  the 
sign  of  7)j jj.  However,  this  term  is  small  for 

and  7’tf1|>r/]|.  Equations  (52) — (54)  indicate  that  both  77/±>0 
and  V^>0  favor  the  excitation  of  the  drift  waves,  whereas, 
77t  >0  and  77^1  >  0  oppose  the  excitation  For  instance,  if  the 
small  quantities,  bn  |(fc1/&n)(  V^0/flr)|,  and  the  ion  Landau 
resonance  term,  are  neglected  in  Eqs.  (52)-(54),  then  we 
have 


Re  (x), 


,  1 2  2  ’ 

1  +  kxps 


(55) 


Ima  /jr  Re  /  k\(f,  ~  1  kVd  \ 
Re  w,  "V  2  kVTr]  y  ]  +  k^  fi  2  ^  «»,  / 


(56) 


Solution  (56)  indicates  instability  (Im6O>0)  for  V^>0  and 
77d|<0.  When  rje{]  >0  and  the  first  term  within  the  parenthe¬ 
ses  of  Eq.  (56)  can  be  neglected  in  comparison  to  77^  /2,  the 


instability  condition  (Im6o>0)  for  the  current-driven  drift 
wave  is 


V^>—  Vel 


(57) 


On  the  other  hand,  for  Vd=Q  (no  current)  and  77,  >0,  the 
instability  condition  (Im6o>0)  is 


1 

-Ve\\< 


i+*2±  p: 


(58) 


If  rfe  =0,  the  instability  condition  is  satisfied  with  any  k2  p2 
>0.  This  is  a  well-known  result  and  the  instability  is  often 
referred  to  as  the  “universal”  instability.  Next,  we  consider 
Uy^kfs^oy^.  Assuming  further  that  k~  p~<gi  1 ,  the  disper¬ 
sion  relation  (50)  simplifies  into 


1  - 


7'iii  *  Vn 


(59) 


where  A  is  the  velocity  shear  parameter  defined  by 


nm)(v;a) 

1  +(*!/*„)(  VW 


(60) 


If  A2>0  so  that  A  is  real,  the  approximate  solution  of  Eq. 
(59)  is 


It  describes  current-driven  ion-acoustic  waves  (in  the  refer¬ 
ence  frame  moving  with  the  ion  flow  velocity)  in  the  pres¬ 
ence  of  velocity  shear,  which  becomes  unstable  (Im6o>0) 
when  Vd  (i.e.,  current)  exceeds  a  threshold  value  Vjj,  where 


(62) 


The  second  term  on  the  right-hand  side,  which  represents  ion 
Landau  damping,  can  be  made  smaller  by  increasing  the 
value  of  and  A,  so  that  the  instability  can  be  excited 

with  a  lower  value  of  V^.  In  the  absence  of  velocity  shear 
(i.e.,  A2=l)  and  for  electron-proton  plasma,  Eq.  (62)  yields 
V°d/c,s  10.13,2.38,  1.17,  and  1.0  for  TefiITin=  10,  15,20,  and 
30,  respectively.  In  the  presence  of  velocity  shear  (i.e., 
A2^  1),  minimizing  Eq.  (62)  with  respect  to  A  we  find  that 
the  minimum  value  of  V^/c5  is  obtained  for  A  =  Am,  where 
Am  is  to  be  determined  from 


1  + 


3- 


exp| 


_Zk 

27, 


=  0. 


(63) 
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For  electron-ion  plasma,  Eq  (63)  yields  Am=  1.478,  1.221, 
1.066,  and  0.88  for  Te  /  7^=10,  15,  20,  and  30,  respectively, 
i.e.,  Am  increases  as  T^ITtl  decreases.  Consequently,  ion 
Landau  damping  for  smaller  values  of  T^IT^  is  reduced  by 
the  velocity  shear.  Substituting  these  values  in  Eq.  (62)  we 
find  the  minimum  values  of  V^/c5=  1.56,  1.28,  1.12,  and 
0.92  for  7^=10,  15,  20,  and  30,  respectively.  Since 
\(k1/kl)(VfeQ/Sle)\  is  typically  negligible,  A  can  be  made 
larger  than  unity  for  |V^/fi(|  <  1  by  adjusting  the  values  of 
kL/k .  Comparison  with  the  homogeneous  (no  velocity 
shear)  case  indicates  that  ion  velocity  shear  allows  excitation 
of  ion  acoustic  instability  with  significantly  smaller  currents 
for  lower  values  of  T^IT^  and  that  the  threshold  current 
does  not  have  as  strong  a  dependence  on  T^IT^  as  it  is  in  the 
homogeneous  case.  This  benign  effect  of  velocity  shear  was 
first  pointed  out  by  Gavrishchaka  et  al .,20  who  showed  by 
numerical  solution  of  the  dispersion  relation  (26)  (without 
the  density  and  temperature  gradients)  that  the  threshold  cur¬ 
rent  for  ion-acoustic  instability  is  significantly  smaller  in  the 
presence  of  shear  and  that  the  threshold  value  is  almost  in¬ 
sensitive  to  variations  of  T^IT,^  over  a  wide  range  of  values 
of  T^iTjw  (from  0.1  to  10). 

If  A2<0  so  that  A  is  pure  imaginary,  and  if  the  small 
imaginary  terms  are  neglected,  Eq.  (59)  describes  a  purely 
growing  mode  with 

(0i  =  i\A\kuCs.  (64) 

This  is  referred  to  as  the  shear-driven  ion-acoustic 
instability  1 

Another  instability  that  has  been  discussed  in  the  litera¬ 
ture  is  the  ion  temperature  gradient  driven  instability. 

This  instability,  like  the  shear-driven  instability,  does  not  rely 
on  or  is  not  affected  by  Landau  resonance  (wave- particle 
resonance).  It  is  realized  when  77,n^>  1  and  in  the  limits  when 
(50)  is  reduced  to 

k2c2 

1  +  -^w*/%  =  0.  (65) 

w; 

It  has  the  solution 


5,  =  -(l  ±  z\  3) 


,  2  2 
*ii c. 


1/3 


w*i  Vi  * 


(66) 


where  the  upper  sign  (+)  is  for  c >0  and  the  lower  sign 
(-)  is  for  co+jTjt  <0. 

In  the  remainder  of  this  section,  we  examine  the  same 
low  frequency  waves  when  different  kappa  distributions  are 
used  to  describe  the  equilibrium  plasma  state. 

(b)  Kappa-Maxwellian:  The  series  and  asymptotic  ex¬ 
pansions  of  Z^M(s),  which  appears  in  Eq.  (31)  and  which  is 
the  same  as  Z*M(s)  introduced  by  Hellberg  and  Mace,  are 


kuzV(k-  1/2) 

2k-  1  (  2  k+  1  «, 

- m  i  — <  + 

K  \  3  K 


kl  <  i , 

(67) 


/VtjTX/c) 


k{i2T(k-  1/2)  (1 


—  1  + 


2k- 3  sf 


_+  ... 


h|  »  1,  (68) 


for  integer  values  of  k.  For  noninteger  (excluding  half- 
integers)  values  of  k,  the  series  expansion  (67)  remains  un¬ 
changed;  but  the  first  term  in  the  asymptotic  expansion  (68) 
is  modified  as  /— ►/-tan  kit.  For  half-integer  values  of  *, 
one  can  first  relate  Z*M  to  Z*  of  Summers  and  Thome1  [see 
Eq.  (60)  of  Ref.  9]  and  then  use  the  series  and  asymptotic 
expansions  of  Z*  for  half-integer  values  of  k  given  by  Sum¬ 
mers  et  al.  The  result  is  that  the  series  expansion  Eq.  (67) 
remains  unchanged;  but  the  first  term  in  the  asymptotic  ex¬ 
pansion  Eq.  (68)  is  modified  as  /7r— ►  iir+  \og(a/s2),  where  a 
is  some  number.  Thus,  noninteger  (including  half-integer) 
values  of  k  will  add  very  small  corrections  to  the  real  part  of 
the  frequencies  of  the  waves  considered  here.  We  neglect  the 
modifications  and  use  Eqs.  (67)  and  (68)  in  the  following 
analysis.  Using  r0(fc,)«l,  V0(b,)  =  ]-biJ  k\De  <£  1 

(quasineutrality  condition),  and  the  expansions  (67)  and  (68) 
in  Eq.  (31),  and  then  keeping  the  leading  terms,  we  find 


J 


2k-  1  /  k±V[0\  ,  , 

- - -(  1  +  —~( +  *1  Ps  -  — [1  -  *.( 1  +  Vo)]  - 

2k-  $  \  k{]  1 /  co, 
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1  7T-— (1  +  V,\ \-biVi 

k  il,  co, 
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+  i\  - F(k )-^~ 
V  2  ’k  VTe 


»  +  -T7f  1-^ 
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2k- 3  \  k  O., 


'*</  1 


W  1 


+ '  \/  ~F(k)I  “ 
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LYk,  -  . 


_  2k-  3  \  k{ i  fli  /  co. 


k2c;  l 
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co]  2k  -3  co, 
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-(*+1) 


=  0, 


(69) 
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where 


-ML  l-ZO-Y2 

/c1  2T(k-  1/2)  \2/c-3/ 


(70) 


The  other  notations  are  the  same  as  before.  It  may  be  verified 
that  F(k)~+  1  and  Eq.  (69)  — >  Eq.  (50)  as  /c— >oo. 

In  the  drift  wave  approximation  (co^  co*f^>k  cs),  the 
dispersion  relation  becomes  simplified  as 


2k-  if.  k  Y^o\  2  2 


IT  .  .  0)e 

K  V‘ 


Tel I 


Ww  1 

co 


2«  L  ,  LYk 

2k- 3  V  ne 


7 j  *e 

+  i  \  j  F{  k) - 

'  2  %V„ 
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1  +  1  x-^jVm-biVu 


1  5 

2k  -  3  Vjn 


X  1  + 


-2  \-(«+l) 


2k  3A:‘V/^.(1| 


=  0 


(71) 


and  its  approximate  solution  is 


Re  co i  =  B{kL,k  )[l  -  6,(1  +  J7i±)]w*f, 


(72) 


2/c-  3  /:||  Qf  /  V7711 
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1  +1  K  “Mi 


1  /  Re  co^  2 


1  + 


1  I  Re  a>;  ^ 2 


2/C-3U'  V’ 


H  v  7YII  f  J 


2/c-  3  V^iiV^i 
-(*+!)] 


(73) 


B(k^k) 


_ 1 _ 

[(2k-  l)/(2/c-3)][l  +(kjkmne)]  +  k\py 

(74) 


Both  the  real  and  the  imaginary  parts  of  co  are  modified  from 
their  values  for  the  bi-Maxwellian  distribution.  The  ion  Lan¬ 
dau  damping  term  in  lm  co  (third  term  within  the  curly 
bracket)  has  a  power  law,  instead  of  exponential  dependence 
on  Re  tbf/knVjin  when  Re  coil  k\{VTi{^>  1  and,  hence,  is  larger 
than  that  for  the  bi-Maxwellian  distribution  in  the  parameter 
regime  Re  c otlk  VY,-^  1.  For  the  limiting  conditions  that  led 
to  Eqs.  (55)  and  (56)  we  now  find 


Re  co,- 


(2 k  -  1  )/(2 k  -  3)  +  k\ p2  ’ 


(75) 


Im  co 


Re  co, 


Re  co,- 
^11^7'dl 


*2Pj2-1/(2k-3) 
(2k-  1)/(2k-  3)  +  p2 


I  2k  kVd 

—  n,s  + - 

2/e  2k -3  co*f 


(76) 


A  comparison  of  Eq.  (75)  with  Eq.  (55)  shows  that  Re  c is 
reduced  from  its  value  for  the  bi-Maxwellian  distnbution,  as 
(2k-  1 ) / (2/c— 3)  >  1.  When  77^.  >0  and  the  first  term  within 
the  parentheses  of  Eq.  (76)  can  be  neglected  in  comparison 
to  7jf  /  2,  the  instability  condition  (Im  co>0)  for  the  current- 
driven  drift  wave  is 


Vd> 


(2k- 3 
\  2/c 


CO 


*e 


(77) 


A  comparison  with  Eq.  (57)  shows  that  the  threshold  current 
for  the  instability  is  reduced  from  its  value  for  the  bi- 
Maxwellian  distribution  by  the  factor  (2/c- 3) /2/c.  In  the  ab¬ 
sence  of  current  (V^=0)  and  for  rjt  >0,  the  instability  con¬ 
dition  is 


]_  k\p]-  l/(2/c-3) 

2 77,11  <  (2k-  1  )/(2k  —  3)  +  k\p2s 


(78) 


If  =0,  the  instability  condition  is  satisfied  for  k 2  p2 
>l/(2/c-3),  which  is  more  stringent  than  the  condition 
k2p]> 0  obtained  for  the  bi-Maxwellian  distnbution. 

Next,  we  consider  the  ion-acoustic  waves  for  which 
co*,.  Assuming  further  that  fc2p2<^Cl,  the  disper¬ 
sion  relation  that  follows  from  Eq.  (69)  is 


2/c-  1 
2/c-  3 


-A2 


,k2c: 


■  +  * 


2/c 
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„  k\\VTe 


•  + A2 


co, 


7/ii  k  Vm 


1  + 


2/c-3A2V2.; 


-(*+0 


=  0 


(79) 


and  its  approximate  solution  for  A2>0  is 


Re  co,  =  KJc  cs. 


(80) 


lm  co 
Re  col 


-(«+!) 


(81) 


where 


a;= 


/2/c- 3 
V  2/c  -  1 


A2. 


(82) 


We  again  note  that  both  the  real  and  the  imaginary  pans  of  co 
are  modified  from  their  values  for  the  bi-Maxwellian  distri- 
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bution.  In  particular,  the  power-law,  instead  of  exponential, 
dependence  of  the  ion  Landau  damping  term  (second  term 
within  braces  in  Im  co)  on  7^/7, when  1 ,  results  in 

larger  ion  Landau  damping  and  changes  the  instability  con¬ 
dition  significantly.  The  threshold  value  for  instability  is 
now  given  by 


Cs 


K  + 


2  k-  1 
2k -3 


1  e\\  A  2 

—A. 


3/2 


X 


1  + 


A2. 

2k-3tJ 


(83) 


In  the  absence  of  velocity  shear  (A2=l),  A2=(2 k 
-3)! {2k-  1)  <  1.  If  we  take  k=  3,  for  example,  then 
A2  =  3/5.  For  electron-proton  plasma,  Eq.  (83)  yields  V^/cs 
a  13.73,  8.31,  5.52,  and  3.05  for  7,]|/7,11=  10,  15,  20,  and  30, 
respectively.  For  a  smaller  (larger)  value  of  k ,  the  values  of 
V°dlcs  are  larger  (smaller).  A  comparison  with  the  corre¬ 
sponding  results  for  the  bi-Maxwellian  distribution  shows 
that  the  threshold  currents  are  larger  for  the  kappa- 
Maxwellian  distribution  due  to  the  larger  ion  Landau  damp¬ 
ing  rates. 

In  the  presence  of  velocity  shear  (A2^  1),  minimizing 
Eq.  (83)  with  respect  to  AK  we  find  that  the  minimum  value 
of  V^/cs  is  obtained  for  AK=AKm ,  where  A^  is  determined 
from 


L  (  AL  TjAk+2  2/c-l/m,\1/2/7-,|\3/2 
\  2k-3TJ  2K-3\mJ  \  Tj 
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I  1hz1LiA2 

\  2k  -3  r,n  " 


-3 


(84) 


Using  Eq.  (84)  in  Eq.  (83)  the  minimum  value  of  V^/cs  is 
obtained  as 


VS)  [2/c<2/c-3)](r  T,)A;m-2  , 

cj ™n  [(.2k -\)/(2k-  3)](7’f[/7’lll)A^m  -  3  *m' 


For  electron- proton  plasma  and  for  k=  3,  for  example, 
Eq.  (84)  yields  AKm=  1.986,  1.681,  1.496,  and  1.269  for 
7^||/7j-j|=  10,  15,  20,  and  30,  respectively,  i.e.,  AKm  increases 
as  7,,/7,-y  decreases.  Substituting  these  values  into  Eq.  (85) 
we  find  the  minimum  values  of  V^/c5  =  2.43,  2.06,  1.83,  and 
1.55  for  7^/7,11=10,  15,  20,  and  30,  respectively.  For  a 
smaller  (larger)  value  of  at,  the  minimum  values  of  V°d/cs  are 
larger  (smaller).  Thus,  in  the  presence  of  velocity  shear,  the 
threshold  currents  for  the  ion-acoustic  instability  are  reduced 
more  significantly  for  smaller  values  of  7e!l/7,1;,  as  in  the  case 
of  the  bi-Maxwellian  distribution,  but  the  reduced  values  are 
still  larger  than  those  for  the  bi-Max wellian  distribution.  This 
again  is  due  to  the  increased  ion  Landau  damping  rates  re¬ 
sulting  from  the  kappa-Maxwellian  distribution. 

The  previously  mentioned  shear-driven  ion-acoustic  in¬ 
stability,  which  is  excited  when  A2<0  (i.e.,  A2<0),  is  now 
described  by 

"<  =  i|Aj*i|CJ.  (86) 

When  compared  with  the  bi-Maxwellian  result  [see  Eq. 
(64)],  the  growth  rate  of  the  instability  for  kappa-Maxwellian 


distribution  is  reduced  by  the  factor  [(2ac—  3) / (2ac—  1)]1/2  for 
the  same  value  of  |A|. 

Referring  to  Eq.  (69),  the  ion  temperature  gradient 
driven  instability,  mentioned  above,  is  described  by 


2k-  1 
2k -3 


+ 


,  2  2 
k  ct 


iVi\  i 


0 


(87) 


with  the  solution 


/  \  3) 


2/c-3  tccj 

2k  ~  1  (w*,%)2 


1/3 


(88) 


The  upper  sign  (+)  is  for  >0  and  the  lower  sign  (-)  is 
for  a>*, 77,11  <0.  When  compared  with  the  bi-Maxwellian 
result  [see  Eq.  (66)],  growth  rate  of  the  instability  for 
kappa-Maxwellian  distribution  is  reduced  by  the  factor 
[(2k-3)I(2k-  l)]1  \  with  other  parameters  remaining  the 
same. 

(c)  Product  Bi-Lorentzian  (PBL):  The  dispersion  func¬ 
tion  ZfL(s),  which  appears  in  Eq.  (38),  is  similar  to  Z*(q)  of 
Summers  and  Thome.  The  series  and  asymptotic  expansions 
of  Z*(s)  for  integer  and  half-integer  values  of  k  are  given  in 
Refs.  1  and  25,  respectively.  Mace  and  Hellberg  extended 
their  results  to  noninteger  (excluding  half-integer)  values  of 
k  The  series  expansion  of  Z*(s)  is  same  for  all  real  values  of 
k,  and,  as  explained  above  [Sec.  IV  (b)],  the  modification  of 
the  asymptotic  expansion  of  Z*(s)  due  to  noninteger  (includ¬ 
ing  half-integer)  values  of  k  may  be  neglected  in  our  analy¬ 
sis.  We,  therefore,  use  the  expansions  for  integer  values  of  k\ 
given  by 


Z?L(sf)  = 


/  \  7rF(  AC  +  1)  /  K  +  1  2 


— Ji 

*!l  V 


2k  +3  \ 


K\  < 1 . 

(89) 


2fL(s,)  = 


Z\  7tT(aC  +1)  1 

*|/2r(/f  +  i/2)(i+s?/#fB)*+1 


1  + 


Kj  \ 

2k  -  1  s2 


|s/|  »  1 .  (90) 


Since  be^  0,  we  keep  only  the  p= 0  term  in  the  electron  sum 
in  Eq.  (38),  and  in  order  to  retain  effects  of  the  order  of  b ,•  we 
keep  p= 0  and  p-  1  terms  for  ions.  Using  £o.o=U  8o.i=~2 
[see  Eq.  (34)],  k\De  <  1,  and  keeping  the  leading  terms  we 
find 
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Here, 


n<c  +  n  (  2*1,  V'2 

“  '  /fa1/2r(K„+  l/2)\2#c1(-  1  / 


(92) 


and  the  other  notations  are  the  same  as  before. 

The  simplified  dispersion  relation  in  the  drift  wave  ap¬ 
proximation  (<0;—  is 


2k  +  l/  k1V,e0\  7  2 

:( '  +  T7T l+k-LPs~—  n-^0  +  %i)] 

2K\\  1  \  A^||  /  (Of 


2Ki  +  2(l  +  k1Vt0 


+  ‘\  .  .  . 

2  2/f„  -  1  ^ii  ii, 

(O 


“>*'  1 
1  -r%ii 


2  *  Vrj|| 


1+(^ii+2'77'1'-^77'i 


j 

-K  -  1 


x  1  + 


J<£ 

2/c  -  1 


-(*,1+2) 


=  0 


(93) 


and  its  approximate  solution  is 


Re  (Dj  =  C(*j.,A:||)[l  -b,(\  + 


(94) 


7 T  Re  CD;  I  /  1 

Imws  |(l--77d jw* 

2k+2/  ^v;0\  _  Vra 

1  +  —  ~  Re  oje  -  <o*e - 


2  k  -  1  \  k  Cle 


1  +  (k  +  “  |  Vi\\  ~  hi  Vt 


i/ 

vm 

l  i  Re  (Oj  \ 2 


1  +  ■ 


2  k  -  1  U,VTln 


1  /  Re  (Of 


2k  -  \  \  k\ym/  j 
(95) 


[(2k  +  l)/(2 k  -  1)][1  +  (k1/kl)(V'e0/ile)]  + k2±p2s 

(96) 

Both  the  real  and  the  imaginary  parts  of  co  are  different  from 
those  for  the  bi-Maxwellian  distribution  as  well  as  for  the 
kappa-Maxwellian  distribution  The  power-law  dependence 
of  the  ion  Landau  damping  term  (third  term  within  the  curly 
bracket  in  Im  a>)  on  Re  (d,fk  Vm ,  when  Re  (otlk  vra»l,  is 
also  different  from  that  for  the  kappa-Maxwellian  distribu¬ 
tion.  This  leads  to  an  instability  condition  that  is  different  not 
only  from  the  condition  for  the  bi-Maxwellian  distribution 
but  also  from  the  condition  for  the  kappa-Maxwellian  distri¬ 
bution.  As  in  the  cases  of  bi-Maxwellian  and  kappa- 
Maxwellian  distributions,  we  ignore  |(A:1/&1|)(V^0/fle)|  and 
the  effects  of  finite  ion  temperature  (i.e„  we  neglect  b(  and 
ion  Landau  damping  term).  Then, 


Re  (Of 


(2k  +  1  )/(2/c,|  —  \)  +  k\p]' 


(97) 


Im  co 


Re  (Of 


k]Ps~  l/(2«  -  1) 

(2 *!+  l)/(2*r  -  1)  +  k\p] 


1  2k  +2  k  Vd 

~77,  + - 

2  /e  2k  -  1 


(98) 


A  comparison  of  Eq.  (97)  with  Eqs.  (55)  and  (75)  shows  that 
Re  a),  is  reduced  from  its  value  for  the  bi-Maxwellian  distri¬ 
bution,  as  (2kh+  l)/(2/C||- 1)>  1,  and  the  reduced  value  is 
different  from  that  for  the  kappa-Maxwellian  distribution. 
When  77^11  >  0  and  the  first  term  within  the  parentheses  of  Eq. 
(98)  can  be  neglected  in  comparison  to  77^/ 2,  the  instability 
condition  (Im  cu>0)  for  the  current-driven  drift  wave  is 


V*> 


/  2k  -  1 
\  2k[{  +  2 


(O 


*c 


(99) 


Thus  the  threshold  current  for  instability  is  reduced  by  the 
factor  (2k  -  l)/(2/c11+2),  when  compared  with  that  for  the 
bi-Maxwellian  distribution  [see  Eq.  (57)].  Also,  the  reduction 
factor  is  different  from  that  for  the  kappa-Maxwellian  distri- 
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bution  [see  Eq.  (77)].  In  the  absence  of  current  (V^=0)  and 
for  77,11  >0,  the  instability  condition  is 


\_  ,  *2jn;-l/(2*-l) 

2  T)e  (2*,  +  1  )/(2/f||  -  \  )+k\p] 


(100) 


If  77,u=0,  the  instability  condition  is  satisfied  for  k\p2 
>  1  /(2/Cj|-  1),  which  is  more  stringent  than  the  condition  for 
the  bi-Maxwellian  distribution  and  is  different  from  that  for 
the  kappa-Maxwellian  distribution. 

Next,  we  consider  the  ion-acoustic  waves  for  which 
Assuming  further  that  k2  p2<g.  1,  the  disper¬ 
sion  relation  that  follows  from  Eq.  (91)  is 

2*u  +  1 


2k  -  1 


a2  pL/  2k,\  +  2 

-a  ~r  +  '\  rn^ii); - 


2*1-1 


L^il^u 


X 


X  1  + 


■  +  A 


2  '  HI 


7/ii  ^ii^r/ii 


1 


2*h-1 


-(*11+2) 


=  0 


and  its  approximate  solution  for  A2 >0  is 
Re  w,  =  A  Kkfs, 


(101) 


(102) 


Re  a), 


Imoj,  _  /ir„,  ,2kh  +  2 

- s  a/-F(*„) - 


8 

2k  +  1 


2#Cn  -  1  \  T, 


where 


2*ii  +  1 


—A2 


3/2 


-  (“-A. 

m ,  \  c* 


1  +  • 


2*  -  1  T iii 


-(*11+2) 


(103) 


A?  = 


2*n  ~  1 
2*  +  1 


A2. 


(104) 


Both  the  real  and  the  imaginary  parts  of  co  are  different  from 
their  values  for  the  bi-Maxwellian  distribution  as  well  as  the 
kappa-Maxwellian  distribution  The  power-law  dependence 
of  the  ion  Landau  damping  term  (second  term  within  the 
braces  in  Im  (o )  on  7V,  /  7^,  when  Te  1 ,  results  in  larger 

ion  Landau  damping  rate  than  that  for  the  bi-Maxwellian 
distribution,  and  it  is  different  from  the  power-law  depen¬ 
dence  for  the  kappa-Maxwellian  distribution.  The  threshold 
value  for  instability  is 


X 


2*,+  1  /  ni; 
2/Cj|  —  1  \  mf 


1  +  ■ 


Teli 


2k  -  1  T, 


A*, +2) 


3/2 


(105) 


In  the  absence  of  velocity  shear  (A2=l),  A2=(2 k 
-  l)/(2/c11+ 1).  If  we  take  *h  =  3,  for  example,  then  A?  =5/7. 
For  electron-proton  plasma,  Eq.  (105)  yields  vj/cf=14.4, 
7.7,  4.64,  and  2.29  for  T^IT^-  10,  15,  20,  and  30,  respec¬ 
tively.  For  a  smaller  (larger)  value  of  #cu,  the  values  of  V^/cs 
are  larger  (smaller). 


In  the  presence  of  velocity  shear  (A2^  1),  minimizing 
Eq.  (105)  with  respect  to  A*,  we  find  that  the  minimum  value 
Of  V°/c,  is  obtained  for  A^=A^,  where  \K  fn  is  determined 
from 


2#c„-  1  7,,,/  2k  -  1  \me/  \  / / 
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(— ^AL-3 

\2*|  -  1  Tn  “i" 


A 


2 
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(106) 


With  the  use  of  Eq.  (106)  in  Eq.  (105)  the  minimum  value  of 
V°dlcs  is  obtained  as 

/v§\  [(2*ii  +  2)/(2*  -  1  )](7~gM/7~f|l) Ay^,  —  2 

UL~ [(2*u  +  1  )/(2k  -  l)](7’V7’,,)A2|m-3'V“m' 

(107) 


For  electron-proton  plasma  and  for  *u=3,  Eq.  (106)  yields 
A^s  1.916,  1.614,  1.429,  and  1.203  for  10,  15,20, 

and  30,  respectively,  i.e.,  A**,  increases  as  T^fT^  decreases. 
Substituting  these  values  into  Eq.  (107)  we  find  the  mini¬ 
mum  values  of  V^/cs  =  2.25,  1.89,  1.67,  and  1.40  for 
r,||/7'/1|=  10,  15,  20,  and  30,  respectively.  For  a  smaller 
(larger)  value  of  /cj),  the  minimum  values  of  V^/cs  are  larger 
(smaller). 

A  comparison  of  the  results  for  the  product  bi-Lorentzian 
distribution  with  the  corresponding  results  for  the  bi- 
Maxwellian  and  the  kappa-Maxwellian  distributions  show 
that  the  threshold  currents  for  ion-acoustic  instability  in  the 
absence  of  velocity  shear  are  larger  than  those  for  the  bi- 
Maxwellian  distribution,  but  are  somewhat  smaller  than 
those  for  the  kappa-Maxwellian  distribution.  In  the  presence 
of  velocity  shear,  the  threshold  currents  are  more  signifi¬ 
cantly  reduced  for  smaller  values  of  T^/Tt ,  as  in  the  cases  of 
bi-Maxwellian  and  kappa-Maxwellian  distributions.  The  re¬ 
duced  threshold  currents  are  still  larger  than  those  for  the 
bi-Maxwellian  distribution,  but  are  somewhat  smaller  than 
those  for  the  kappa-Maxwellian  distribution. 

The  shear-driven  ion-acoustic  instability,  which  is  ex¬ 
cited  when  A2  <0  (i.e.,  A2<0),  is  now  described  by 

‘3<  =  '|A*J*i;Cj.  (108) 


When  compared  with  the  bi-Maxwellian  result  [see  Eq. 
(64)],  the  growth  rate  of  the  instability  is  reduced  by  the 
factor  [(2/c;)- \)/(2k  + 1)]1/2  for  the  same  value  of  |A|. 
However,  the  growth  rate  is  somewhat  larger  than  that 
for  the  kappa-Maxwellian  distribution  [see  Eq.  (86)],  as 

|aK||>|aJ. 

Referring  to  Eq.  (91),  the  ion  temperature  gradient 
driven  instability  is  now  described  by 


2*  +  1 
2k  -  1 


t2c2 


(109) 


with  the  solution 
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w;=  -(1  ±  i\'3) 


2k  -  1 
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2k  +  1  (w*,»7fli)2 


1/3 


(110) 
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i\  7rr(/c) 


/c1 2r(/c- 1/2) 
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2k-  1 

• )  - - s. 


The  upper  sign  (+)  is  for  £u*/*7iii>0  and  the  lower  sign  (-)  is 
for  £o*(?7l1  <0.  The  growth  rate  of  the  instability  is  reduced 
by  the  factor  [(2ku- 1)/(2k,|+ 1)]1/3,  when  compared  with 
that  for  the  bi-Maxwellian  distribution  [see  Eq.  (66)],  and  it 
is  slightly  larger  than  that  for  the  kappa-Maxwellian  distri¬ 
bution  [see  Eq.  (88)],  other  parameters  remaining  the  same. 

(d)  Bi-Lorentzian  Plasma:  Since  be~ 0,  we  keep  only  the 
p= 0  term  in  the  electron  sum  in  Eq.  (47).  Consequently,  we 
need  to  know  the  series  expansion  of  Z^q(s^)  for  1.  For 
the  ions,  we  keep  p- 0  and  p-  1  terms  in  Eq.  (47)  in  order  to 
retain  effects  of  order  bt.  Hence,  we  need  to  know  the 
asymptotic  expansions  of  Z^(s,)  and  Z^ffc/)  for  1.  As 
we  mentioned  earlier,  Z®  q  is  same  as  Z*M(s)  and  so  its  series 
and  asymptotic  expansions  are  given  by  Eqs.  (67)  and  (68), 
respectively  [see  the  discussion  following  Eqs.  (67)  and 
(68)].  The  asymptotic  expansion  of  Z^s^can  be  derived 
from  the  asymptotic  expansion  of  Z^q(s,)  by  using  the  rela¬ 
tion  (44).  We  thus  have 
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(113) 


Using  the  expansions  in  Eq.  (47),  assuming  k\De (  <$C  1 ,  and 
keeping  only  the  leading  terms,  as  before,  we  find 
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where 
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The  simplified  dispersion  relation  in  the  drift  wave  approximation  (a>;~ (o^^kpC^  is 


(115) 
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(116) 


We  notice  that,  except  for  the  contribution  proportional  to 
birjiL  in  the  ion  Landau  damping  term,  Eq.  (1 16)  is  identical 
to  the  corresponding  dispersion  relation  for  the  kappa- 
Maxwellian  distribution  [see  Eq.  (71)].  So,  in  the  limit  when 
biVii  ,  the  stability  properties  of  the  drift  waves  for  bi- 


Lorentzian  distributions  are  same  as  those  for  the  kappa- 
Maxwellian  distribution  [see  Eqs.  (72)-(78)]. 

For  the  ion-acoustic  waves  (oj^k  c^o)^)  we  assume, 
as  before,  that  /c^p?<Cl,  and  the  dispersion  relation  that  fol¬ 
lows  is 
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which  is  identical  to  the  dispersion  relation  (79)  for  the 
kappa-Maxwellian  distribution.  Hence,  the  analysis  of  the 
ion-acoustic  wave  dispersion  relation  for  bi-Lorentzian  dis¬ 
tribution  is  same  as  that  for  the  kappa-Maxwellian  distribu¬ 
tion  [see  Eqs.  (80)— (85)].  We  mention  here  that  the  lon- 
acoustic  instability  for  Lorentzian  distnbution  (i.e., 
Ta\\=Tal)  and  in  the  absence  of  velocity  shear  has  been  nu¬ 
merically  analyzed  by  Meng  et  ai  The  analyses  for  the 
shear-driven  instability  and  the  ion  temperature-gradient 
driven  instability  for  bi-Lorentzian  distnbution  also  remain 
the  same  as  those  for  the  kappa-Maxwellian  distribution  [see 
Eqs.  (86)— (88)]. 


V.  SUMMARY  AND  DISCUSSIONS 

We  have  presented  the  linear  dispersion  relations  for 
electrostatic  waves  in  spatially  inhomogeneous,  current- 
carrying  anisotropic  plasma,  where  the  equilibnum  particle 
velocity  distnbutions  are  modeled  by  vanous  Lorentzian 
(kappa)  distribution  functions  and  by  the  well-known  bi- 
Maxwellian  distribution.  Spatial  inhomogeneities,  assumed 
to  be  weak,  include  density  gradients,  temperature  gradients, 
and  gradients  (shear)  in  the  parallel  (to  the  ambient  magnetic 
field)  flow  velocities  associated  with  the  current.  In  order  to 
illustrate  the  distinguishing  features  of  the  kappa  distribu¬ 
tions,  stability  properties  of  the  low  frequency  (lower  than 
ion  cyclotron  frequency)  and  long  perpendicular  wavelength 
(longer  than  ion  gyroradius)  modes  have  been  studied  in  de¬ 
tail,  and  the  results  have  been  contrasted  with  those  for  the 
bi-Maxwellian  distribution.  Specific  attention  has  been  given 
to  the  drift  waves,  the  current-driven  ion-acoustic  waves  in 
the  presence  of  velocity  shear,  the  velocity  shear-driven  ion- 
acoustic  modes,  and  the  ion  temperature  gradient  driven 
modes. 

The  growth  rates  of  the  drift  wave  instability  and  the 
current-driven  ion-acoustic  instability,  both  of  which  rely  on 
wave-particle  interactions  for  their  excitation,  are  reduced 
from  their  values  for  the  bi-Maxwellian  distribution  due  to 
larger  ion  Landau  damping  rates  associated  with  the  kappa 
distributions.  For  the  same  reason,  the  excitation  conditions 
for  these  two  instabilities  are  more  stringent  in  the  case  of 
the  kappa  distributions.  The  dominant  ion  Landau  damping 
rates  are  proportional  to  ImZ'  (prime  on  Z  denotes  deriva¬ 
tive  with  respect  to  its  argument),  which  have  power-law 
dependence  on  Re  <5,/ (A: n  Vn\)  »  1  for  the  kappa  distributions. 
This  is  in  sharp  contrast  with  the  exponential  dependence  for 
the  bi-Maxwellian  distribution.  As  a  result,  the  ion  Landau 
damping  rates  of  plasma  waves  that  are  excited  in  the 
Re  <oi/(fc!|Vyil|);»  1  regime  are  larger  in  kappa-distnbution 
plasmas  than  in  bi-Maxwellian  plasma.  Figure  2  shows  the 
marked  differences  in  the  behavior  of  Im  ZMs  associated  with 
the  different  velocity  distributions,  when  Re  £,/(£  VTi )»  1. 
A  particularly  important  consequence  of  the  enhanced  ion 


Landau  damping  rates  is  that  the  threshold  currents  for  the 
ion-acoustic  instability  in  kappa-distnbution  plasmas  are 
larger  than  those  in  the  bi-Maxwellian  plasma,  even  in  the 
presence  of  shear  in  the  parallel  flow  velocity.  Relativistic 
effects  associated  with  the  suprathermal  ions  that  participate 
in  the  considered  instabilities  have  been  neglected  in  the 
present  nonrelativistic  treatment  under  the  assumption 
(7,||i»7’,1)<Cm/c2.  For  a  recent  paper  on  the  modeling  of 
energetic  particles  by  relativistic  kappa  distribution,  see 
Ref.  27. 

The  stability  characteristics  of  the  other  two  instabilities 
(shear-dnven  ion-acoustic  instability  and  ion  temperature- 
gradient  driven  instability),  which  do  not  rely  on  wave- 
particle  interactions  and  for  which  Landau  damping/growth 
terms  in  the  dispersion  relation  may  be  neglected,  can  be 
better  understood  in  terms  of  (ne\/n0)  and  (niX/n0)y  where  neX 
and  nn  are  the  density  perturbations,  since  the  dispersion 
relations  are  obtained  by  demanding  quasineutrality  (nfX 
=  nn).  (nfl/n0)  and  (/T, j //i0)  are  also  helpful  in  understanding 
the  origin  of  the  reduced  frequencies  of  the  drift  waves  and 
the  ion-acoustic  waves  in  kappa-distribution  plasmas.  Under 
the  conditions  assumed  in  Sec.  IV,  the  electron  density  per¬ 
turbations  for  the  different  equilibrium  distributions  are 


(ne\ln0)  =  a(kL,k ) 
r 

for  BM 

[(2k-  l)/(2K-3)](e#,/rJ,forKM 
[(2k  +  l)/(2*,|-  l)](e^,/rj,  for  PBL 
(Same  as  the  expression  for  KM),  for  BL, 

(118) 


where  a(k^  )=  1  +(A'1/A  |)(V^0/n^).  The  differences  in  the 
expressions  can  be  traced  to  the  differences  in  the  series 
expansions  of  Z  for  <oel(k  Vr^)<C  1  and  to  the  differences  in 
the  relations  of  6^  to  Te .  The  ion  density  perturbation,  on  the 


-ImZ  • 


FIG.  2.  Comparison  of  -ImZ'  vs  at,  where  *  =  Re  io(f(hym).  Solid  curve 
represents  -ImZ8^,  dashed  curve  represents  -lmZ*KM  for  k= 3,  dashed- 
dotted  curve  represents  -Im  Z^PBL  for  k  =  3,  and  the  dotted  curve  represents 
-ImZ^BL  for  *= 3.  As  noted  in  the  text,  Z'k0*l=Z’kkm. 
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other  hand,  is  same  for  all  the  distributions  and  is  given  by 


.  V'-o  \ 

I-T7T-— u  +  %) 

k  it: 


U): 


2„2 


Ti 


<0? 


(O: 


e<t>  i 


(119) 


when  b(  and  k\p?s  are  neglected  for  simplicity.  This  is  due  to 
the  fact  that  the  leading  term  in  the  asymptotic  expansions  of 
Z  for  a5,/(/C||V7>/]|) » 1  is  the  same  for  all  the  distributions. 
Equation  (118)  shows  that  the  adiabatic  response  of  the  elec¬ 
trons  to  the  electrostatic  potential  perturbation  in  kappa- 
distribution  plasmas  is  reduced  from  its  value  in  bi- 
Max  wellian  plasma.  The  reduced  adiabatic  response  of  the 
electrons  explains  the  reduction  of  the  growth  rates  of  the 
shear-driven  ion-acoustic  instability  and  the  ion  temperature- 
gradient  driven  instability  as  well  as  the  reduction  of  the 
frequencies  of  the  drift  waves  and  the  ion-acoustic  waves  in 
kappa-distribution  plasmas. 
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